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Abstract 

A light stop, with an i?-parity- violating coupling X'l^i, has been suggested 
as an explanation of the excess in high-Q^ neutral current events observed 
at the HERA collider. We show that in this scheme a corresponding excess 
in charged current events — such as that reported by the HI Collaboration 
— can appear naturally, without calling for the presence of light sleptons or 
additional i?-parity-violating couplings, if there exists a chargino lighter than 
the stop. The predicted event shapes agree well with the data. The relevant 
region of parameter space is identified, taking into account constraints coming 
from precision electroweak measurements, atomic parity violation and recent 
searches for first-generation leptoquarks at the Tevatron collider. 



CERN-TH/97-176 
July 1997 



The reported [^,0] excess of neutral-current (NC) events at large in deep inelastic 
scattering (DIS) at the HERA collider has been interpreted as a possible hint of physics 
beyond the Standard Model (SM) [@-[l2|. Among the various solutions considered, the 
currently favoured ones involve the s-channel production of a scalar 'leptoquark', which 
couples to the positron and one of the valence (or sea) quarks in the proton. Scalars 
of this nature, namely squarks, form an essential ingredient of i?-parity violating ) 
supersymmetric theories. Supersymmetry (SUSY) not only provides one of the most 
well- motivated extensions of the SM but also leads to an explanation of the observed NC 
excess P-[7|,p!^. Its role as a framework for solutions to the HERA anomaly thus merits 



serious consideration. 



Within the framework of the minimal supersymmetric extension of the Standard 
Model (MSSM) the relevant contribution to the superpotential reads 

W = K,,L,Q,Dl , (1) 

where i, j and k are generation indices, Li and Qj are the SU (2) l doublet lepton and quark 
superfields respectively, while is a charge-conjugate right-handed down quark super- 
field. In view of the various constraints on couplings from low-energy processes ||13| , 



the most natural interpretation of the HERA NC data is given by the s-channel produc- 
tion |]14[ of a left-handed charm (cl) or top (ti) squark, with a mass of (9(200) GeV and 
a coupling to the valence d-quark, 

A'yi ~ = Br(g, ^ e+rf) , (2) 

VPed 

where j = 2, 3. It is perhaps worth mentioning that a sizeable effect proceeding from the 
interaction of the positron with the sea quarks in the proton would, in general, demand 
rather large couplings, which are generally inconsistent with constraints from low-energy 
data. The only exception is the possibility of stop production from a sea s-quark 
through A'i32 — 0A/^/J3^. The parameter space in this scenario has subsequently been 
shown to be constrained by LEP measurements as well as by limits on the electron- 
neutrino mass [|16[ . In the following, we shall concentrate on a valence quark interpretation 
of the HERA anomaly. 

Searches at the Tevatron collider for generic leptoquarks decaying into lepton and 
jet — of which squarks of ^p supersymmetry are one example — put strong restrictions on 
the allowed branching ratio /Jed- For a squark mass of (9(200) GeV, the bounds from the 



CDF and DO analyses can be combined to indicate 

/3ed < 0.5 . (3) 

While this result seems to rule out most explanations of the HERA anomaly involving 
leptoquark resonances^, which have no other decay channels, supersymmetric models 
can easily get around the impasse. Given the rich particle spectrum of SUSY models, 
R-paritj- conserving interactions can induce additional decay processes (not counting, of 

^Attempts have been made, however, to construct models that evade this constraint p9[|. 
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course, those that lead to hard charged leptons plus jets in the final state). If this is indeed 
the case, with the negative results at the Tevatron being due to an i?-parity-conserving 
decay that is unobservable with the current data sample [^, then it is natural to ask 
what is predicted at HERA, where a major fraction (1 — Ped) of the produced squarks will 
decay through these modes. One should typically expect some additional signals beyond 
those expected within the SM. 



As a matter of fact, the HI Collaboration [|T| — and recently, the ZEUS Collabora- 



tion as well — have reported a small excess in charged current (CC) events at large 
in the same data sample that contains the NC excess. Although less significant sta- 
tistically than the latter, the observed number is nevertheless too large to be consistent 
with SM background calculations at 2.5 standard deviations. In view of the expectation 
that i?-parity-conserving decays of the produced squark would be observed at HERA, an 
identification of the CC excess with these modes seems to be the most natural conclusion. 
It is thus interesting to ask whether SUSY scenarios can admit this possibility without 
contradicting other known phenomenological results. Such an investigation assumes par- 
ticular importance since it is difficult to accommodate the CC excess within other 
non-supersymmetric explanations for the NC events. Some attempts in this direction 
have already appeared in the literature 0,|T2|,^,|23|. In Ref. [Q, it was pointed out that 
i?-parity- conserving decays of the squark could possibly lead to the CC signal. However, 
the observed event shape — low jet multiplicity, large missing momentum, relatively large 
hadronic invariant mass of the final state, and, of course, a high — need to be taken 
into account in any serious pursuit of such an explanation. A more detailed study of the 
event shape has, in fact, been undertaken, mainly in the context of charm-squark pro- 



duction, in Ref. |T^, where an explanation of the CC excess is obtained by postulating 
a light slepton-sneutrino pair, each of which decays into jets through couplings. In 
order to obtain a signal consistent with the data, especially with the low jet multiplicity^, 
the sneutrinos must be rather light, mj> < 80 GeV, while the wino mass M2 should not 
be much larger than 200 GeV. Similar requirements lead to a CC signal consistent with 
the experimental data in the context of top-squark production |1T2|,[2^. Although this 
is a plausible explanation of the HERA anomaly, light sneutrinos and low values of M2 
lead to a sizeable sneutrino pair-production cross section at LEP2 pSf and consequently 
to observable rates for four-jet final states This last prediction presents a knotty 



problem, in view of the difference between results presented by the ALEPH Collaboration 
and the other LEP experiments but the general hope is that a definitive result may 
be obtained from the next run of LEP2. It is also worth pointing out that the solution 
suggested in Ref. |]T2[ may lead, for scharm production, to events with a small fraction 



of wrong-sign i~+ jet final states at HERA, which could be detectable ||28[ as more data 
accumulate. A final state with fi~^+ jets (with or without missing momentum) can also 
be accommodated in their scenario in the presence of a light smuon-sneutrino pair. 



In this letter, we describe a scenario in which a CC excess with the appropriate event 
shape may be generated even in the absence of light sleptons. We concentrate on the stop 

^Although, even then, the CC events would, in general, tend to have two visible jets (see Fig. 2 of 
Ref. IH). 
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interpretation with an Ifip coupling \'^^^ consistent with Eq. (j^). A hght stop in the 200 
GeV mass range is natural in models with radiative breaking of electroweak symmetry 
where squarks of the first two generations are heavy ^9) . It is perhaps better motivated, 
therefore, than the competing option of a light charm-squark. The hierarchy of masses of 
the left-handed stop and the sbottom arising from the SU (2)-breaking condition, 

m?^ — m?^ = ml — ml + cos 2/5 > , (4) 

suggests a simple way of generating CC events: if an i?-parity-conserving decay mode 
of the stop leads to sbottoms in the final state, the CC events may be generated by the 
subsequent decay vd oi the sbottom through the same I/ip coupling X'l^i- (A similar 
scenario for the cl is not possible, since the sl would be somewhat heavier than the cl 
for tan [3 > 1.) 



In principle, if the sbottom is light enough (m^ ^ 120 GeV), a substantial branching 
fraction for the weak decay of a stop into a bottom squark and a (with the last going 
mainly into jets) may be allowed. This has been proposed recently by Kon et al. as the 
cause of the CC excess at HERA. Although t ^ h + —>■ jets is a straightforward 
and economical scenario, it has several potential problems. First, it would mostly lead to 
CC events with a jet multiplicity of two or three, whereas, of the four events observed by 
HI, only one is consistent with more than a single jet [00]. Moreover, the hadronic mass 



distribution is rather broad and one has to appeal to the low statistics p3[ to claim any 
agreement with the data. Apart from such kinematic arguments, one also encounters the 
indubitable fact that the contribution of the stop-sbottom doublet to the p parameter 
tends to be too large (Ap^ > 2 x 10~^) if the stop-sbottom mass splitting is so large. 
Slightly smaller values of Ap may be achieved by introducing large mixing angles in the 
stop and sbottom sectors, and tuning them appropriately. This has, in fact, been invoked 
by the authors of Ref. [E^]. However, even then, Ap^ > 1.5 x 10~^ and hence the scenario 
is only marginally consistent at the 2a level [E4|. In addition, constraints from atomic 



parity- violation experiments tend to be more severe in the case of large stop mixing |31 



All these difficulties can, of course, be overcome trivially by pushing the sbottom mass 
to larger values. This would immediately make the jets coming from the particle 
softer, and hence less often visible, while increasing the missing pt and the invariant 
hadronic mass of the CC events. It also reduces the contribution to the p parameter. 
Unfortunately, when the sbottom mass is pushed high enough to reduce the contribution to 
the p parameter to acceptable values, the i?-parity-conserving branching fraction becomes 
very much smaller than the one, and the corresponding CC event rate gets highly 
suppressed. 



It is thus clear that a phenomenologically consistent explanation for the CC obser- 
vation will require light supersymmetric particles in addition to the {t, b) pair. The light 
slepton option has been commented on above and will not be looked into any further. Since 
the stop decay into neutralinos is naturally suppressed^, the obvious course is to have a 

^One could consider a non-negligible mixing between stop and the left-handed charm squark, and thus 
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chargino lighter than the stop. For most of the MSSM parameter space, this also means a 
neutralino that is considerably lighter and is, therefore, the lightest supersymmetric parti- 
cle (LSP). Assuming that there are no other light sfermions, just three decay channels are 
now open to the chargino, namely into the neutralino {x~^ x^W* — > x'^ff), into the 
sbottom (x"*" cb) through the Cabibbo-Kobayashi-Maskawa (CKM) mixing (provided 
the sbottom is lighter than the chargino), and an decay through a virtual stop. For 
most of the parameter space of interest, however, the decay of the chargino has very 
small width. We shall assume, in the following, that the t-b mass splitting is < 35 GeV, 
which seems optimum for the p-parameter constraint, and that > > m^. Now, 
depending on the branching ratios, two distinct possibilities present themselves. 

(i) The chargino decays predominantly into the neutralino. If the latter is lighter than 
the sbottom (m^o < mi), it can have only a three-body decay — >■ bdu^, bdue)- In this 
case, the final state would typically have multijets and relatively low missing momentum. 
Such configurations will probably be lost in the DIS (CC) background coming from SM 
interactions. On the other hand, if the neutralino is heavier than the sbottom (which 
thus becomes the LSP), i.e. 

mi > > m^M > m-^ , 

then it would decay into b + bi and the Ijip two-body (and only) decay of this b^ would 
result in a significant amount of missing momentum. Although, at first sight, it might 
seem that we would still have too many jets, most of these are now very soft as long as 
mi — m^ < 35 GeV. A simple parton-level Monte Carlo simulation shows that most of 
the jets (apart from the one coming from the bi decay) fail to meet the minimum energy 
criteria of the HI experiment for defining a jet [0,0 — this typically results in most events 
ending up with only a single 'jet' and almost none with more than two. Even when the 
parton energy is adequate, the second jet would be soft and broad, and it is debatable 
whether it would be detectable at all. It is quite easy to confirm that the kinematic profiles 
for these configurations match those for the observed CC events. Thus, a cascade decay 
of stop into chargino into neutralino into sbottom seems to be an attractive scenario. 
However, like the ones considered above, this too suffers from an inherent drawback. 
Since all the charginos produced in the stop decay (and we must remember that this is 
the dominant decay of the stop for I3ed ~ 0.5) are potential CC event candidates (modulo 
detection efficiency), we now tend to get too many of the latter. Unless one wishes to 
appeal to experimental effects and the low statistics, this number can be reduced only 
by providing the sbottom with another decay channel. What can that be? If we violate 
gaugino mass unification and choose the soft supersymmetry-breaking parameters such 
that the mass of the lightest neutralino (x?) is well below that of the sbottom while the 
next-to-lightest neutralino (^2) lies just below the chargino and is the one involved in the 
above chain, then, indeed, a new decay mode 6l — >■ 6x1 can be opened. A scan of the 
parameter space shows that this rather complicated scenario can be achieved for a limited 
part of the parameter space and that too, only if there are large deviations from gaugino 

allow the tree-level decay ti, — > cx°. It can be shown, however, that the mixing angle must be very large 
(~ 1) to lead to the observed signal; such a large mixing angle is difficult to accommodate within any 
reasonable theoretical framework. 
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mass universality. Though this scheme works — as we have checked — it has hmited 
aesthetic appeal and will not concern us any further. 

{ii) The chargino decays predominantly into the sbottom and a charm quark. This 
mode is naturally suppressed by the smallness of the CKM matrix element Vcb = 0.036 — 
0.046. Although some enhancement can be obtained in the event of a nonzero b-s mixing, 
we prefer to be conservative and do not consider this option. Thus for —>■ cbi to be 
the dominant decay mode, the chargino-neutralino coupling which drives the competing 
decay (x+ f f) has to be quite small. Fortunately, this occurs naturally when 
is large compared to M2. The lightest neutralino is then mostly bino-like and couples 
very weakly to the PF, while the lighter chargino and the next-to-lightest neutralino are 
mostly wino-like and nearly degenerate in mass. Thus, the chargino decay into the second 
neutralino will be strongly suppressed, if not disallowed. The mass hierarchy obtained in 
this scheme is of the form 

mi > ~ m^o > mi> m^o . 

We find that a branching ratio (3ed ^ 0.5 can naturally be obtained for values of M2 ~ 
~ 190 GeV) and corresponding values of X[^-^ as determined by Eq. (||). In order 
that the jet resulting from the charm quark should not be observable, we also require the 
sbottom to be fairly heavy, typically in the 170-180 GeV range, which is already indicated 
by the p parameter constraint. The decay of the sbottom into ud now provides a hard 
jet and missing momentum, as in the cases discussed above, which can be the origin of 
the observed CC excess. 

For this last scenario — which will concern us in the rest of this letter — we observe 
that the sbottom can also decay in the i?-parity-conserving mode — >■ bxi, with the 
neutralino then decaying into 

Xi bdue , bdue . (5) 

through an intermediate sbottom. If the neutralino is light enough, m^-^ ~ 100-140 GeV, 
the missing of this decay mode will tend to be low, and will mostly fail the CC selection 
criteria imposed by the HI Collaboration (see below). Potential CC events resulting from 
this decay channel will then be lost in the SM DIS background. We can thus have a 
reduced number of CC events without making any assumption beyond that of a large \fj,\ 
and M2 ~ 190 GeV. The simplicity and economy of this scheme is apparent. In effect, we 
only assume: 

A. a small stop-sbottom mass-splitting, which is easily achieved with a modicum of 
left-right mixing in the stop sector; 

B. a chargino, which is mainly wino-like, lying between the stop and sbottom in mass; 
and 

C. a light neutralino, which is dominantly bino-like, and whose mass is considerably 
below that of the sbottom, typically 100-140 GeV. 
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With these assumptions, a study of the predicted event shapes and efficiencies is now in 
order. 



In the presence of left-right mixing = imphes that the hghter eigenstate is a 
pure ii), the stop production cross section at HERA is given by 

= YdfEo ^^'^31 -ose.f K u (^) , (6) 



where E^, ES, are, respectively, the initial electron and proton energies in the laboratory 



frame and K ~ 1.3 parametrizes the QCD correction I^Sj. For fdix), the c?-quark distri- 
bution within the proton, we use the MRS(H) structure functions, calculated using 
the package pdflib In this analysis, we neglect effects due to initial state radiation. 



The NC and CC signals are estimated using a parton-level Monte Carlo event generator. 
For jets within the angular coverage of their detector, the HI Collaboration requires a 
minimum transverse momentum of 

f^T>l5 GeV . (7) 

In our parton-level analysis, each final-state quark is assumed to give rise to jets. All such 
'jets' that lie within a cone 



ARjj = ^At]"^ + A02 < 1.0 (8) 

of each other are 'merged' to form a single jet by the simple expedient of adding their 
four-momenta together vectorially and identifying the sum with the momentum of the 
combination. Here Ar/ is the difference of pseudorapidities and Acj) denotes azimuthal 
separation. 

For NC events, the HI analysis requires that the final-state positron must lie within 

10° <e^< 145° , (9) 

(the proton direction is positive) and it must have a minimum transverse energy: 

PTe > 25 GeV . (10) 

In addition, the missing transverse momentum must satisfy ]1)t / a/Pte < 3 GeV^/^ The 
DIS Lorentz invariants can be estimated using either of the pairs (^e,^jet) and {9e,Ee); 
the latter give : 

l/e = l-§sin2^e, Ql = -^, M,= 

1 - Ve 

where E^ = 27.5 GeV is the energy of the incoming positron beam. The NC events are 
required to satisfy 

0.1<i/e<0.9. (11) 
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With the above set of cuts (and other standard ones required to remove cosmic-ray back- 
grounds) the HI Collaboration reports 12 events with Ql > 15000 GeV^ for an inte- 
grated luminosity of 14.2 pb~^ against the expected DIS background of 4.7 ±0.76 events. 



For the CC events, on the other hand, the cuts @-(|TT]) are no longer operative. 
Instead, the HI analysis demands that the missing transverse momentum be large: 

i)T = Y. PT(hadrons) > 50 GeV . (12) 



Analogous to the NC case, one can define DIS variables ||3^ using only measurable 
hadronic variables: 

nE-p\\) ^2 i>T 



2i?o ' l-Vk' 

where the sum runs over all observed hadronic clusters (in our simulation, 'jets' that pass 
the cuts). With a selection criterion of 

Vh < 0.9 , (13) 

the HI Collaboration observes 4 events above Ql > 15000 GeV^, whereas the expected 
SM background is 1.77 ± 0.87 events. 



It is easy to see that the efficiency for the NC events should be solely determined by 
the stop mass m^-, and for ~ 200 GeV, can be approximated by 

€NC ~ 0.56 , (14) 

where we assume detector efficiency to be unity. The CC efficiencies, on the other hand, 
can depend on all the relevant masses in the decay chain. Again, our computations show 
that the only significant dependence is on the bottom-squark mass, and 

ecc ~ 0.37 (0.45) ± 0.008, for = 170 (180) GeV . (15) 



As the NC event topology has been studied elsewhere |1^ , and shown to be 
consistent with the data, we shall not repeat the exercise here. In Fig. we present the 
kinematic distributions for the CC events. To be concrete, we have chosen to depict it 
for = 200 GeV, = 190 GeV and for three typical sbottom masses. It turns out 
that the phase-space distribution is not very sensitive to the exact value of my.+ and is 
primarily governed by the sbottom mass. 

Each of Fig. |l](a-c) depict the px distribution of the three potential 'jets' for a given 
mi. As expected, the hardest jet emanates mainly from the 6-decay, but with (small) 
contributions from the soft jets that merged with it. Consequently, its p^-distribution 
just reflects the Jacobian peak with a slight smearing due to the small transverse boosts. 
For the next-to-hardest jet (dashed lines), the situation is quite different. In general, 
it has relatively low px and is mostly rejected by the cut (|^). However, depending on 
the t-x^ si-iid the x^-b mass splittings, some conflgurations may have sufficient pt to be 
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Figure 1: Event topology for mi = 200 GeV, = 190 GeV, and for different values of 
mi- (a-c) The pt distributions for the three 'jets' in each case, (d) The Q\ and (e) the 
Mh distributions for all three cases. 



observable. The remaining jet, if it has not merged with the hardest two, is always too 
soft to be seen. This is amply demonstrated by the dotted lines. Thus, even though we 
started with as many as three quarks in the final state, we never have more than two jets 
at the observable level. Moreover, the fraction of 2-jet events is always much smaller than 
that of single-jet events for ~ 180 GeV; the fraction naturally goes up as the sbottom 
mass decreases. For the case where two jets can indeed be seen, their angular separation 
tends to be evenly distributed within the interval 1 < ARjj < 3. 



Figure ^{d) shows that the number of CC events is a rather sensitive function of Ql; 
most of the events are concentrated within ISOOOGeV^ < Q\ < 25000GeV^. This is well in 



agreement with the present HERA data Events with Ql > 30000 GeV^ 



2111 would, 



however, be harder to explain within this scenario. The mass distribution in Fig. |l](e) 
peaks at about 180-190 GeV, but has a long tail, which may account for the fact that 
one of the observed events has a low Mh (~ 157 GeV). It is interesting that such a peak 
in the mass distribution can be obtained in the scenario of Ref. only for rather small 
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sneutrino masses, while the scenario of Ref. ||2^ has a much flatter and broader mass 
distribution. 



In Fig. |I] we present results for only a few points in the available parameter space, 
but it is easy to see that the same qualitative features would hold for other points as long 
as these satisfy the conditions (A-C) above. What of the events resulting from the decay 
chain of Eq.(|^)? A simulation analogous to the one described above shows that most of 
these events indeed fail to survive the cuts of Eqs.(^ |T^-|TB|), with the missing momentum 
cut proving to be the most severe. The corresponding efficiency is a sensitive function 
of m^o and ranges from about 4% to 13% as the latter varies from 90 GeV to 140 GeV. 
This is rather too low to yield an observable number of multijet events with large missing 
momentum in the present data sample. 



Having established that the HERA results for CC events can be explained within the 
present scenario, it becomes interesting to identify the part of the MSSM parameter space 
that supports our solution. Certain requirements are immediately obvious. As explained 
above, since we want the chargino to be wino-like and slightly lighter than the stop, we 
are immediately constrained to M2 ~ 190 GeV. For n, tan/3 and the ratio M1/M2, the 
situation is more complicated. We note that the ratio of the number of NC events to CC 
events in this scenario is given by 



with efficiencies as in Eqs. (14) and ( p3|) . Consider a typical value f3ed — 0.4 consistent 
with the Tevatron constraints. Since Br(6 vd) k, 0.4 (0.6) for m^o ~ 100 (140) GeV, we 
have Ncc/Nnc — (0.35-0.5) xBr(x+ ^ cb). A good agreement with the HI rates thus 
requires Br(x+ cb) > 0.5. This constrains the parameter space in the /x-tan/? plane, 
and we present the favoured region in Fig. |^. While this region depends on the value 
of the coupling X[^^ as well as on the squark masses (these not only determine the 
production cross section and the efficiency, but also bound M2), we present the parameter 
space for only one representative set of values. To be deemed acceptable, a point in the 
parameter space was required to lead to 7 ± 1 NC events and 2.5 ± 1 CC events for an 
integrated luminosity of 14.2 pb~^ and simultaneously satisfy 0.3 < Ped < 0.5. Note that 
these bounds on P^d are relaxed for a slightly heavier stop (say 210 GeV). This would 
considerably expand the favoured region. A similar effect can be obtained by considering 
a mixing the b-s sector. Our estimate should thus be regarded as a conservative one. 



As shown by the solid lines in Fig. ^ for our scenario to work within the scheme of 
universal gaugino masses (M1/M2 = | tan^ — 0.5), a relatively large is indicated. 
The curious shape of the curves owes its origin to the existence of minima in the wino 
content of the LSP as tan (5 is varied. The asymmetry between positive and negative fi 
can be traced to a similar source. While /i > 1 TeV (or /x < —600 GeV) is perfectly 
acceptable, it is interesting to ask if smaller values could be accommodated. An elegant 
way to achieve this is to reduce the mass difference —m^o so that the two-body decay 
X+ X^W~^ is suppressed kinematically and the relevant process is the three-body decay 
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Figure 2: The favoured range (see text) in the /i-tan/3 plane for {X'l-^i = 0.07, m^- = 
200 GeV, = 175 GeV) for four values 0/M1/M2. We assume a small stop mixing. 



of the chargino into the LSP and two hght fermions. This happens naturally in the event 
of small corrections to the universality relation in the gaugino sector. In Fig. ^ we also 
show the favoured region for three other values of M1/M2. A priori, the dependence on 
this ratio could be quite complicated. While a larger value of M1/M2 would kinematically 
suppress the ~^ x'^ff channels, it would also increase the higgsino content of both 
the chargino and the LSP and consequently the x~^x^W coupling. At the same time, the 
X'^cb coupling decreases too. Still, for the range of interest, the dependence is monotonic. 
It must be noted though that M1/M2 > 0.7 implies m^o > 140 GeV, which runs the risk 
of resulting in a few (possibly observable) multijet events with large missing momentum. 



In Fig. ^ we have chosen particular values for the squark mixing. A non-zero value 
for 61 will serve to modify both Br(x"'" cb) and Br (6 vd). Since the bbB coupling 
vanishes exactly for tan^^ = 0.5, for a mixing angle in the vicinity of this value, the 
curves in Fig. |^ would move inwards by approximately 150 GeV. Stop mixing, on the 
other hand, manifests itself mainly in determining the stop-chargino coupling and hence 
I3ed- For \9i\ < 15°, the dependence on this parameter is negligible for negative /i. For 
positive the deviation from Fig. |] is maximum for M1/M2 ~ 0.5 and, for moderate 
9l, could, at best, lead to a shift by ~ 100-150 GeV, with the weakest constraints being 
obtained for 6^ ^ 0. 



A further interesting feature of this scenario is the possibility of the chargino decaying 
into a neutralino and a i'^u pair instead of jets. This would typically lead to a hard lepton, 
jets and missing momentum in the final state. For i = fi, this could perhaps explain ||36[ 
the muonic event observed by HI. A similar event with the r"*" leads to a narrow 
jet and is swamped by the CC DIS signal. An with low missing momentum would, 
similarly, be indistinguishable from the NC DIS, but events with larger missing momentum 
might be detectable as statistics multiplies. 
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We now turn to the various low-energy measurements that, in principle, could further 
constrain the allowed parameter space. As we have seen in the above discussions, if the 
stop is responsible for both the NC and the CC excesses at HERA, it cannot be a pure 
left-handed state. The corrections to the p parameter would then be too large to be 
consistent with precision electroweak measurement data. While a large mixing can evade 
this problem, such a solution runs counter to the bounds from atomic parity violation |^ 



The effective weak charge measured experimentally now receives a contribution from both 
stops; for a nucleus ^ Xa, it is given by 



AQ 



w 



-{2A-Z) 



1 tan2 6'r 
1 

1 ' 2 



:i7) 



where T denotes the heavier stop. The experimental bound on AQvf of the Cesium 
atom, along with the size of the NC excess at HERA, thus translates into a bound in the 

mi ~ 200 GeV implies, at the 95% C.L. level. 



, rrif) plane. For example. 



I3ed ^ 0.3 



1 + tan^ Oz —k- 



m. 



It should be remembered though that the error in AQ^y is dominated by theoretical 
uncertainties, and so this bound should be considered with great caution. It is clear that, 
given the upper bound on I3ed imposed by the Tevatron data, small values of the mixing 
and/or a large hierarchy between the stop masses are preferred. 



A simultaneous resolution of the p-parameter and the atomic parity-violation con- 
straints thus call for a small t-h mass splitting as well as a small stop mixing angle. While 
this may seem difficult at first, in reality, it can arise in a natural way. Consider the case 
where mu-, the supersymmetry-breaking mass parameter in the right-handed stop sector, 
is large in comparison with that in the left-handed sector (wig). The lightest stop mass 
is then approximately given by 

m?^m| + Di + m? {l-^ , (18) 

where At = At — fi/ tan (3 and = 0.35M| cos2/5. The lighter stop will be mainly 
left handed {i.e. 6^ is small), and the heavier stop mass {ruf ~ mu) may be pushed 
to large values. Clearly, the smaller Oi is, the higher mf is pushed — see Fig. |(a) — and 
consequently, the smaller the additional contribution to ^Qw is. Ignoring the small 
mixing in the sbottom sector, the left-handed sbottom mass will be given by 

m?^m| + D^ + m2, (19) 

where = — 0.42M| cos 2/3. It is clear that, under these conditions, for values of 
~ mu, the lightest stop and sbottom masses will be close and hence the contribution 
to the p parameter will be small. On the other hand, for a fixed positive value of m^ — mi, 
the contribution to the p parameter is smaller for larger values of mu and for ~ mu- 
This conclusion remains valid as long as the mixing in the sbottom sector is small. Large 
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Figure 3: (a) The heavier stop mass as a function of the mixing angle for = 200 
GeV, mi = 175 GeV, Ifi] = 1 TeV, A,, = 0, tud = 2 TeV and mu < 1 TeV. (b) The 
corresponding contribution of the stop-sbottom doublet to the parameter Ap. 

9i tends to push the value of mq up, leading to a substantial mixing in the stop sector for 
the same values of the stop masses. Sbottom mixing appears naturally in the large tan j3 
region and its effects are clearly seen in Fig. |5(a). 



In Fig. 0(&), we plot the stop contribution [^] to Ap as a function of 6i for 



mi = 200 

GeV and = 175 GeV. Fixing fi and mo, the right-handed sbottom mass parameter, 
allows us to see the effects of a small mixing in the sbottom sector for large values of tan /3. 
Large values of fi and tan (3 also have a striking effect on the branching ratio Br(6 —>■ sj) 
as we shall discuss below. As argued above, the larger the left-handed component of the 
lighter stop, the smaller Ap. The experimental numbers suggest Ap^ < 1.5 x 10^^ 
for mt = 175 GeV and mh — 100 GeV. From the figure, it is apparent that 6i ^ 15° is 
preferred by the data. Observe that even lower values of Ap may be obtained by relaxing 
the upper bound on m^j. 



Since both stop masses are determined in this scenario, it is necessary to discuss 
their impact on the mass of the lightest CP-even Higgs [Q. In Fig.^(a), we present 
the variation of the Higgs mass as a function of the mixing angle, for the same values 
of the stop and sbottom masses. Although there are four branches associated with the 
two different signs of p and the mixing parameter At, for > 250 GeV, there is little 
difference between the branches. Hence, we plot the Higgs mass only for positive values 
of p and At. We see that small mixing angles are preferred for low values of tan/3, while 
no information on the mixing angle may be obtained from Higgs mass considerations in 
the large tan f3 regime. 



Since the stop and chargino masses are fixed, we can determine the stop-chargino 
contribution to Br (6 57) It is important to remember though that the physical 
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Figure 4: (a) The mass of the lightest CP- even Higgs as a function of the stop mixing 
angle, for the same parameters as in Fig.^ and = 1 TeV. (b) Brffe s'y) as a 
function of the stop mixing angle for the same parameters. The branches are defined by 
(sgn(/i),sgn(A)). 



branching ratio depends on additional parameters which are not relevant for an inter- 
pretation of the HERA data. For a specific choice of these parameters and m£+ = 190 
GeV, Fig. §(6), shows the branching ratio as a function of the stop mixing angle. We set 
the first and second generation squark masses to 2 TeV and consider all four branches 
defined by the signs of fi and At. The QCD scale has been fixed at Q = mb/2, a value 
that reproduces the next-to leading-order correction^ in the SM. Setting Q = mh would 
typically reduce the rates by about 20%, giving a measure of the theoretical uncertain- 
ties in these calculations. For tan/5 > 10, an acceptable rate may be obtained only by 
large cancellations between different contributions to the amplitude, which demands spe- 
cific values of the Higgs and supersymmetric mass parameters. For the above parameter 
choice and = — ITeV, we find no solutions for tan/3 = 10, 50. Similarly, for /i = +lTeV 
and tan/? = 50, acceptable solutions may be found for only one of the branches of A^. 
Although this seems to prefer lower values of tan/3, the situation can change for other 
parameter choices. In fact, with suitable parameters, we can always reproduce the exper- 
imental results for Br (6 57) in the scenario under consideration. 

In all of the above, we have held mn = 2 TeV. Although our scenario is consistent with 
lower values of this mass parameter, the analysis becomes more involved in the large tan (3 
region, where the mixing in the sbottom sector becomes large. As noted earlier, a large 
value of 9i will modify the relevant branching fractions, allowing for lower values of /i for 
a given M1/M2 and tan/?. On the other hand, since large 9i generally implies a large 6*^-, 
the atomic parity violation constraints become more stringent in this case. Observe that, 
for an appropriate choice of the mass parameters, larger values of the sbottom mixing 

complete next-to leading-order computation in the supersymmetric model has not yet been per- 
formed. 
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angle may allow for a reduction of the value of Ap, even for moderately large values of 
the stop mixing angle. 

To summarize, then, we have investigated the possibility that both the neutral and 
the charged current anomalies seen at HERA at high are consequences of the resonant 
production of a stop, which decays to an e+o? final state to give the NC events and to hy^ 
followed by cascade decays to give the CC events. Something of this nature is indicated 
by the negative results of the Tevatron search for leptoquarks/squarks in the 200 GeV 
mass range. Among the various scenarios discussed, the CC events are best explained 
through the decay cLl, which is CKM-suppresscd, but can be sizeable in the limit 

of a somewhat large \n\ parameter, when the competing decay channel is considerably 
reduced. The kinematic distributions in this scenario explain the observed CC excess 
rather well and the parameter space which supports this solution is perfectly consistent 
with electroweak precision measurements at LEP and with other constraints such as those 
originating from radiative S-dccays. atomic parity violation and Higgs searches at LEP. 
Although it may seem premature to demand precise agreement with data for the CC 
excess reported by HI and ZEUS, our scenario makes firm predictions which should be 
testable in forthcoming runs of the existing high-energy accelerator facilities. 

The authors would like to acknowledge discussions with G. Altarelli, G. Giudice, 
S. Lola and P. Zerwas and arc grateful to Y. Sirois for specific information regarding the 
HI observations. We thank P. Chankowski for confirming our numerical results for Ap^. 
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